Ionization fronts in negative corona discharges. 
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In this paper we use a hydro-dynamic minimal streamer model to study negative corona discharge. 
By reformulating the model in terms of a quantity called shielding factor, we deduce laws for the 
evolution in time of both the radius and the intensity of ionization fronts. We also compute the 
evolution of the front thickness under the conditions for which it diffuses due to the geometry of the 
problem and show its self-similar character. 
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A common feature in transient discharges preceding 
dielectric breakdown is the creation of a non-equilibrium 
plasma through the propagation of a nonlinear ioniza- 
tion wave into a previously non-ionized region. Modern 
concepts of pattern formation which have been already 
applied in different contexts (see e.g Q and 0), have 
also been used in order to gain new analytical insight in 
this old problem 0. 

When a sufficiently high voltage is suddenly applied 
to a medium with low or vanishing conductivity, extend- 
ing fingers of ionized matter develop. They are called 
streamers and are ubiquitous in nature and technology 
A minimal streamer model, consisting of a fluid ap- 
proximation with local field-dependent impact ionization 
reaction in a non-attaching gas like argon or nitrogen, 
has been used to study the basics of streamer dynamics 
00000 EtJ • The essential properties of planar fronts 
for this minimal model have been obtained as a first step 
towards an effective interface description. For the pla- 
nar fronts, the mechanism of dynamical front selection 
has been understood and explained 0, 0. The dis- 
persion relation for transversal Fourier-modes of planar 
negative ionization fronts has been derived in an at- 
tempt to study how the electric screening layer influences 
their stability. 

As a further step, we consider in this paper the evo- 
lution of negative ionization fronts in a gas under the 
influence of a non-uniform external electric field. The 
field is created by a potential difference Vq applied be- 
tween an electrodes pair. The geometry of the electrodes 
determines the non-uniformity of the electric field. Dis- 
charge develops in the high field region near the sharper 
electrode and it spreads out towards the other electrode. 
This type of discharges is called corona. It is a negative 
corona discharge when the electrode with the strongest 
curvature is connected to the negative terminal of the 
power supply. We will consider this case. 

The so-called minimum model consists in the follow- 
ing dimensionless set of equations (for physical param- 
eters and dimensional analysis, we refer to the previous 



discussions in 



^-V-j=n e /(E) ) 

drii 

V • E = m — n e . 



(1) 

(2) 
(3) 



The equation JTJ describes the rate of change of the local 
dimensionless electron density n e . It is equal to the diver- 
gence of the local electron current density j plus a source 
term n e /(E) representing the generation of electron-ion 
pairs due to the impact of accelerated electrons onto neu- 
tral molecules of gas. The value of /(E) is given by the 
Townsend approximation 



/(E) =£exp (-!/£), 



(4) 



where E is the modulus of the local electric field E. In 
equation @ we consider that the rate of change of the 
ion density m is equal to the source term due to impact, 
since we take the ion current density to be negligible in a 
first approximation (the speed of ions is typically much 
smaller than that of electrons). The local value of the 
electron current density is specified as 



j = n e E + DVn e 



(5) 



using Ohm's law in the first term and considering diffu- 
sion effects in the second one. Note that this expression 
does not include the effect of the magnetic field created 
by the motion of electrons, as it is supposed that their 
speed is much smaller than the speed of light. Equation 
(0) is Gauss' law in local form, coupling the electric field 
to the charge densities. 

Since our primary goal in this paper is to address the 
effects of curvature in front propagation, we will neglect 
diffusion effects as in [tj. That allows us to reduce the 
set of equations CD - © into a simpler form in order to 
give analytical results for the evolution of the ionization 
fronts. From Q, an d © with D = 0, we obtain 



d_ 

at 



(m - n e ) + V • (n e E) = 0, 



(6) 



2 



and from (|3J|, taking the time derivative, 

v -{-m) = di {ni - ne) - 

Equations JBJ) and J7J) give then 

V • + n e E] = 0. 



(7) 



(8) 



The term inside brackets in JSJ is, due to Maxwell equa- 
tions, proportional to the curl of the magnetic field in the 
gas. As it is supposed that the magnetic field is negli- 
gible, we can take it equal to zero and integrate in time 
yielding 



E(r,i) = E (r)exp 



/ drn e (r, 
Jo 



r) 



(9) 



which gives the local electric field E in terms of the initial 
electric field Eo and the electron density n e integrated 
in time. Equation @ motivates the definition of the 
quantity 



u(r, t) — exp 



(10) 



If this quantity is completely determined in a particular 
problem, then using equations JHJ, (|10|l and J3J all the 
physical fields can be obtained through the expressions 



E(r,f) 
n e (r,t) 

rii(r,t) 



Eo(r)tt(r,*), 
1 du(r,t) 
u(r,t) dt 

1 du(r,t) 
u(r,t) dt 
V-(E (r)«(r,t)), 



(11) 
(12) 

(13) 



in which the initial condition Eo for the electric field 
should be known. Equation Hll|) reveals clearly the role 
played by the function u(r, t) as a factor modulating the 
electric field E at any time. The electronic density is 
positive so u(r, t) decays damping the electric field. For 
this reason we call it the shielding factor. The shielding 
factor determines a screening length which changes with 
time: a kind of Debye's length which moves with the 
front leaving behind a neutral plasma. 

The problem is thus reduced to finding equations and 
conditions for the shielding factor u(r, t) from equations 
and conditions for the physical quantities E, n e and n,-. 
Substituting l|ll(l - l|13() into the model equations 0-10, 
after some algebraic manipulations and integrating once 
in time, the evolution of u(r, t) is given by 

-— = V • E u - mo(r) - exp — )ds, 14 
u at ./ R „„ \ s I 



I E u 

u(r,0) = u (r) 



1, 



(15) 



where Eo is the modulus of Eo and n^o is the initial ion 
density. Boundary conditions should be imposed depend- 
ing on the particular physical situation. 

In what follows we will consider a typical corona ge- 
ometry: two spherical plates with internal radius Rq and 
Ri » Rq, respectively. An electric potential difference 
Vb is applied to these plates, so that V(R\) — V(Ro) = 
Vb > 0. The initial seed of ionization is taken to be 
neutral so that 



n e o(r) = n i0 (r) = p (r). 



(16) 



We consider the evolution of negative ionization fronts 
towards the positive plate at r = R\. The initial electric 
field Eo(r) between the plates is 



E (r) 



C 



u r , C = V 



RqRi 
Ri — Ro 



(17) 



We substitute l|17|l into equation l|14fl and change the 
spatial variable r to 



3C' 



(18) 



so that the evolution for the screening factor takes the 
form 



On 



du 

u— = -upo(x) - u 
ox 



C V/3 



(^t) 

(g&) 



exp 



-1 



ds. (19) 



The equation l|19|) governing the behaviour of the screen- 
ing factor is a Burgers' type equation, where po(x) is the 
initial distribution of charge. The condition for the ini- 
tial value of the screening factor is, by l|15|) . u(x, 0) = 1. 
Following the usual procedure of resolution of Burgers' 
equation we can integrate equation 1)19(1 along the char- 
acteristics x c (t) defined by 



dx c (t) 
dt 



= u(x c (t),t), 



(20) 



transforming (|19|l into an ordinary differential equation. 

First the case of sufficiently localized initial conditions 
is considered. More specifically, the initial electron den- 
sity strictly vanishes beyond a certain point. Under simi- 
lar conditions, the existence of shock fronts with constant 
velocity has been predicted for the simpler planar geom- 
etry ||. 

Taking a homogeneous thin layer of width S << (i?i — 
Rq) from r = Rq to r — Rq + 5, the initial charge distri- 
bution is then 

n e o(r) = n i0 (r) = p , R < r < R a + 5, 
n e0 (r) = n i0 {r) = 0, R + 6 < r < R x . (21) 

In figure ^ w e show the electron density distribution n e 
which corresponds to some arbitrary choice of parame- 
ters poj S, Ro,Ri,Vq. The electron density has been cal- 
culated using expression i|12|) and plotted as a function of 
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FIG. 1: The shock wave development at regular intervals of 
time when the initial charge distribution is well localized. The 
amplitude of the front calculated analytically is plotted in 
dashed line. 



r at different times t. There appears a sharp shock with 
decaying amplitude, separating the region with charge 
and the region without charge. 

From these numerical data, we can measure the veloc- 
ity of propagation of such front. In figure [21 it is plotted 
the position of the shock rj as a function of time. The 
velocity of propagation is clearly not constant. However, 
if we plot the position of the front in terms of x, one can 
observe the following linear relation (see inset figure 0) 



Xf(t) = t + x a , 



(22) 



which implies, in terms of the original variable r, an 
asymptotic behaviour 



r/(t) ~ (3C) 1/3 t 1/3 



(23) 



for the position of the front. 

Remarkably we can deduce expressions for both the 
amplitude and propagation velocity of the shock in ex- 
plicit analytical form. In order to do that, we write lo- 
cally near the front the solution as 



i(x, t) = 1 — a(t)ip(£). £ = x — Xf(t) 



(24) 



We substitute this ansatz into (|19|) and since the integral 
term is small when x » 1 we get 

a (t)^(0-a(i)^(C)4(i)-a'(i)^(e)+a 2 (t)^(0^(C) « 0, 

(25) 

x f (t) = t + x , (26) 
a(t) = . 13 . , (27) 

<p(£) = p-\x-t-x Q ), (28) 



implying that 



FIG. 2: Position of the shock front versus time. Squares are 
numerical results and the continuous lines are the theoretical 
predictions. In the inset it is shown the position in the scaled 
variable x. It can be observed the linear dependence in this 
variable as predicted in l|26|l . 



where (3 is an arbitrary constant to be fixed from initial 
conditions. Equation lj26|l proves that the position of the 
shock front follows the law l|23[) for spherical geometry. 
From H12(l . the electron density reads 

r i i+( x -t-x )/(t+t ) rr<f , rr 

n e {x,t)&{ *+*° i-(*-t-xo5/(t+to5 ' x ^ l + x o ^ 

I , X > t + Xq 

(29) 

which implies that the amplitude of the front decays as 



n e (x f (t),t) 



1 



t + t Q 



(30) 



In figure ^ the analytical curve (|30|l has been plotted 
in dashed line, showing an excellent agreement with the 
numerical data. 

We want to conclude this paper with a brief discussion 
of the case in which the initial charge distribution is not 
localized as, for instance, one such that 



n e0 (x) = Uio(x) = pa(x) - e" 



x » 1. 



(31) 



For the planar case it was predicted Q a constant velocity 
for the propagation of the front, although no shock front 
would develop unless the decay is sufficiently fast. 

As we did above, we solve the problem numerically 
assuming spherical symmetry, so x is related to the ra- 
dial coordinate r by l|18l) . In figure [3] the solution for 
the electron distribution is shown. The shock front does 
not appear in this case. Instead, a front with increasing 
thickness propagates. In the scaled variable x, we have 
checked that the centre of this front moves with constant 
velocity as the shock front does. These facts are apparent 
from the figure. 
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FIG. 3: Front development when the initial condition is not 
localized, plotted at regular intervals of time. 

Using scaling arguments, it can be shown that the 
asymptotic local behaviour near the front can be de- 
scribed in the following self-similar form: 

"M^)' (32) 

where £ = x — t, f is some universal self-similar pro- 
file and S\ is a constant measuring the front thickness in 
rescaled units. Its value depends on the physical param- 
eters and initial conditions. Hence the front presents a 
typical thickness 

£ c « S x t. (33) 

The fact that, even neglecting diffusion, the front spreads 
out linearly in time is a remarkable feature of the curved 
geometry considered here. For this reason it can be 
termed as geometrical diffusion. In figure^Jwe have plot- 
ted the numerical solutions rescaled according to 
showing a clear convergence towards a universal profile. 

The principal results and contributions from the work 
presented in this paper can be summarised as follows. 
First we have introduced the shielding factor as the fac- 
tor damping the electric field in non-equilibrium elec- 
tric discharges when the magnetic field can be consid- 
ered negligible QlOfl. This factor defines a characteristic 
length analogous to Debye's length for stationary dis- 
charges. The physics contained in the minimum model 
for streamer discharges can be reduced to the study of 
the evolution of the shielding factor. We have derived the 
equation which governs its evolution (|T4"|) for a gas like ni- 
trogen or argon without taking into account the diffusion 
of charged species and the processes of photoionization. 

Then we have consider the case of a negative corona 
discharge with spherical symmetry. In this case, the dis- 
charge takes place in a non-homogeneous electric field 




FIG. 4: Front profiles rescaled according to the self-similar 
law given by expression 1321 . The profiles converge asymp- 
totically to function / in that expression. 

and the equation for the shielding factor turns out to be 
a Burgers' one. We have extended the results of planar 
fronts to this case where the geometry is curved. Depend- 
ing on initial conditions for the charge distribution, one 
might have negative shocks or spreading fronts. In both 
cases, the amplitude decreases in time and the propaga- 
tion velocity follows a power law. In the case of spreading 
fronts we have proved the appearance of diffusion-type 
phenomena due to purely geometrical effects. 
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